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Abstract 

We analyze the thermodynamical potential of a lattice gas model with three 
components and five parameters using the methods of Catastrophe Theory. We 
find the highest singularity, which has codimension five, and establish its transver- 
sality. Hence the corresponding seven-degree Landau potential, the canonical 
form Wigwam or Aq, constitutes the adequate starting point to study the overall 
phase diagram of this model. 
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1 Introduction 



The study of phase diagrams of complex systems is an important part of Thermody- 
namics and applied physics. Phase diagrams can be constructed either from experi- 
mental data or from theoretical models (for example, molecular models). Among these, 
the wide class of lattice gas models is particularly suitable for a mathematical analy- 
sis. Here we will focus our attention on the lattice gas model for a system with three 
components which simulates, in particular, a binary fluid mixture. A wide literature 
has already been devoted to it from different points of view JT], |2], |3|, [|, and references 
therein]. We are mainly interested in attempting to present an overall analysis of its 
phase diagram, with particular attention to its highest multicritical point, that is, the 
one with the highest codimension in the five- dimensional parameter space that we 
consider. 

In the mean field theory, the Gibbs potential is a function of the concentration of 
two of the three components and depends on three thermodynamical parameters, which 
can be taken as the temperature and the chemical potentials of the two components, 
and on three molecular parameters. The phase diagram deduced from this function is 
an accurate description of the system, except close to the (multi) critical points, where 
fluctuations become important and alter significantly the mean field theory predictions. 
For this reason, the Gibbs potential has been the basis for determining the overall 
phase diagram |2|, |3|. The method used in 0] establishes the qualitative features of 
the phase diagram, namely, the instability and (multi) critical (hyper) surfaces which 
divide its various regions, as well as some coexistence (hyper) surfaces. Those methods 
are considerably powerful but not sufficiently rigorous by mathematical standards. 
However, a well established mathematical theory for the analysis of singularities of 
potentials and hence the associated phase diagrams does exist, namely, Catastrophe or 
Singularity Theory (CT). 

Catastrophe theory has indeed been applied to the description of phase transitions 
|| |6| and, in particular, of phase diagrams of complex thermodynamical systems includ- 
ing fluid mixtures || . The philosophy behind these applications is however somewhat 
perverted: one starts with a system on which some knowledge of the phase diagram 
is available, perhaps its salient features, and surmises a polynomial (Landau) poten- 
tial from among the variety supplied by the CT classification (called canonical forms), 
which is supposed to embody the properties of the thermodynamical potential of the 
system near the phase transition of interest. Usually, this potential is well analyzed in 
the mathematical literature and its other properties can be safely assigned to the phys- 
ical system. In summary, this procedure is phenomenological in nature and amounts to 
fitting of phase diagrams. Although it utilizes CT, it is only to take advantage of well 
studied potentials. However, this way is, of course, less powerful and precise than the 
approach determining the polynomial potentials from a complete singularity analysis 
of the actual thermodynamical potential. It is this second approach the one we shall 
adopt here. A big advantage of our approach is that it is presented as a well defined 
algorithm leading to a systematic way to analyze any general problem susceptible to 
be studied within the context of CT. 
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Futhermore, this second alternative agrees with the methods of Ref. 0, which in 
fact approach those of CT from the point of view of classical thermodynamics. (See 0] 
in regard to the convergence of both techniques). We intend to take the best of both 
worlds for the problem in hand: to draw intuition from thermodynamical methods 
and mathematical soundness from the theorems of CT. In particular, we emphasize 
the study of transversality of the actual thermodynamical potentials which guarantees 
that those simple forms (polynomial potentials or canonical forms) represent indeed 
up to a diffeomorphism the original thermodynamical potential. We shall give a brief 
account of the CT algorithm utilized throughout the paper together with an Appendix 
for more mathematical details. 

Our results essentially agree with and support those in ||. However, we hinge 
less on the visualization of the phase diagram and more in the classification of its 
singularities, relying for the construction of the phase diagram on the straightforward 
method of gluing patches, each described by a standard canonical form for which 
the phase diagram can be found in the literature. Besides, we clearly establish the 
possibility of Landau potentials in two variables, that is, of corank-2 canonical forms, for 
the system with three components. This possibility was dismissed in 0. Nevertheless, 
this case is sufficiently complex on its own to postpone it for future work. 

This work is organized as follows: In section 2 we introduce the fundamental con- 
cepts of CT, the more thechnical points of which are left for an appendix. In section 
3 we describe the thermodynamical potential to be analyzed, give its physical inter- 
pretation and discuss general stability questions which help connect usual concepts in 
Thermodynamics with those in Catastrophe Theory. In section 4 we apply the CT pro- 
gram to the potential previously introduced, reducing it first to a one- variable potential. 
In section 5 we show a solution with the highest codimension and its transversality, 
thereby concluding that it is the Wigwam catastrophe. We also study in this section a 
peculiar singularity of lower codimension, associated to the physics of critical azeotropy. 
The last section is devoted to a discussion of the previous results and of the structure 
of the phase diagram entailed by them. 

2 Generalities about Catastrophe Theory 

In this section we are going to review very briefly the main concepts of CT. The reader 
desirous of more technical details is referred to Refs. J7], §] . As is well known, CT deals 
with the singularities of smooth real- valued functions. The nature of these singularities 
is revealed by perturbing those functions. If as a result of a perturbation the qualitative 
properties of the function remain unaffected we will say that this function is stable or 
structurally stable. In other words, a function is said locally stable at a given point if 
there is a smooth change of coordinates so that the new function has the same structure 
as the old function. If a function is stable at all points then we will say that this function 
is globally stable. In a more precise way, a given function is a perturbation of another 
one at a given point if the distance between both functions is arbitrary small. The 
concept of distance leads us to Topology. We can define the Taylor-series topology 
in the space IR d where d gives the number of the Taylor series coefficients. Thus the 
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fc-jet of a given function at a given point is the Taylor series truncated beyond terms 
of degree k. Several definitions of distances can be given and all of these are expressed 
in terms of the fc-jets of each function. 

Now the next important question is what information is lost when we truncate the 
Taylor series of a function around a given point, namely, the problem of determinacy. It 
consists of determining whether a function can be truncated and if so, for what value of 
the degree of the Taylor expansion it can be truncated without any loss of substantial 
information. Furthermore, to determine the most general family of functions of the 
smallest dimension d which contains the original function is called the problem of 
unfolding. The unfolding dimension is the number of parameters describing a general 
perturbation and the minimum number to describe it is called the codimension. When 
all the unfolding terms go to zero, the remainder of the universal unfolding is called 
the germ of the canonical form. 

The next step is to introduce the concept of transversality as a means to study 
structure stability and genericity. This concept was originally introduced by R. Thorn 
|J and, in general, is not widely used to classify physical phenomena in terms of ele- 
mentary catastrophes. A property is called generic if the subset for which the property 
is valid is open and dense in the original set. In other words, when a property is in- 
variant under a perturbation, this property is called generic or structurally stable. The 
theorem of transversality shows that it is a generic property of functions to have only 
isolated critical pointdQ and such functions are stable under perturbations. Two mani- 
folds of H n intersect transversally if either their intersection is empty or they intersect 
transversally at all their points of intersection, that is, the direct sum of their tangent 
spaces at the point has dimension n or they span the tangent space IR n at that point. 

CT has usually not been applied in a rigorous way using all these concepts and 
theorems needed for its correct implementation. We claim in this work that if this is 
done so, the procedure proposed by the theory is not by any means cumbersome and 
time consuming. The Catastrophe Program proposed here provides a very useful and 
systematic way to examine with not very much effort general behaviors of physical 
systems. 

Let F(x, A) be a real function with state variables x±, . . . , x n (x G H n ) and control 
parameters Ai, . . . , A r (A G IR r ); that is, F : IR n+r — > IR. We are to proceed as follows: 

1. We pick (xq,\q) such that Xq is a degenerate critical point of F(x,X) and we 
consider the unfolding f(x, A) = F(x + xo, A + Ao) — F(xo, Ao) and h(x) = f(x, 0). 

2. One calculates the determinacy and codimension of h from the fe-jet of h (see 
Appendix). Of course, if h is fc-determine then h ~ j k {h), that is, the function 
h is equal to j k (h) up to a change of coordinates and hence they are equivalent 

1 We warn the reader not to confuse thermodynamical and mathematical terminology. In thermo- 
dynamical terms a critical point is just an equilibrium point and the word critical is reserved for a 
higher singularity. 
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and have qualitatively the same properties; therefore, cod(/i) = cod(j fc (/i)). 



3. One studies the fc-transversality of F and if this function is fc-tranversal we can 
affirm that F and the canonical form of the unfolding of h are isomorphic and we 
can replace the original F function for this canonical unfolding. If not, we can 
claim that the F function is not susceptible to be studied by CT. 



3 Description of the Gibbs potential 

According to Ref. || a phenomenological model for a ternary mixture is obtained by 
assuming the Gibbs potential in the form 

G = N[a'yz + b'xz + c'xy + R T (x In x + y In y + z In z)\ , (1) 

where N = N x + N y + N z gives the number of total moles and N x , N y and N z the moles 
of each component. The variables x, y, z are the mole fractions defined by x = N x /N, 
y = Ny/N and z = N z /N; and hence we have the constraint 

x + y + z = l (2) 

where < x, y, z < 1. Finally, a', b' and d are phenomenological energy parameters. 
The energy part is the most general quadratic term, given that x+y+z = 1. This model 
can be derived from the mean field theory of a lattice model Hamiltonian with variables 



taking three different states, representing the molecules of the three components ||10 
Then a', b' and d represent molecular interaction parameters. Let us consider the Gibbs 
potential Eq. (JJ) in a reduced form, dividing by NRT, and thus 

G(x, y, z,a,b,c) = ayz + bxz + cxy + xliax + ylny + zlnz (3) 

where now the new parameters a, b, c are defined with respect to the old ones a',b', d 
dividing them by RT. The concentrations are supposed to be determined by some 
boundary conditions, such as the values of the chemical potentials of two compo- 
nents, say fi x and fi y . The mean field theory prescription is then to minimize the 
non- equilibrium Gibbs potential G — \x x x — fi y y with respect to x and y to obtain the 
equilibrium conditions 

dG 8G _ 

~dx~~^~dy~~ ^ 

They allow to solve for x and y as functions of ^ x ,fi y and the parameters a, b and c, 
provided that the Jacobian det ^g7~j^ = det dfjG is not zero. 

Thermo dynamical stability further requires that the matrix dfjG be positive defi- 
nite. This property is called convexity and must hold for any thermodynamical poten- 
tial, except on the instability hyper- surfaces, which are the simplest singularities we 
may encounter in a phase diagram. An instability can occur only near a phase tran- 
sition, when two equilibrium states, one unstable — hence unphysical — and the other 
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meta-stable coalesce and disappear. In other words, a meta-stable state becomes un- 
stable and, consequently, we speak of instability. This is the kind of sudden change 
in the configuration of a system to which CT owes its name. In mathematical terms 
we say that the critical point (equilibrium state) is degenerate. As a consequence, one 
cannot solve for x and y as functions of \i x ,\i y or the solution is multi- valued, corre- 
sponding to the existence of various equilibrium states. The simplest instability occurs 
when only one eigenvalue of the stability matrix vanishes; that is, the instability only 
affects one variable. One is to focus on this variable, which is called relevant, to con- 
sider further singularities. Therefore, it is convenient to transform the potential into a 
function of just the relevant variable by solving the equilibrium conditions for the other 
variables and substituting for them. The standard thermodynamical procedure that 
performs this operation is the Legendre transform. In fact, this procedure can be used 
for potentials in other fields, whenever there is an underlying geometrical structure 
in the total space of variables, including state and control variables, called a contact 
structure. We address the reader interested in the general formulation to the literature 



TTJ. How the Legendre transform is implemented in our case will be seen in the next 
section. Further singularities are studied afterwards with the one-variable potential. 
Next comes what can be called critical instability, followed by the tricritical point and 
so onwards. 

Several systems of interest are described by this Gibbs potential Eq. ([!]): A ternary 
mixture at constant volumen, for example, a mixture of metals; a spin lattice where 
the molecules have spin one; a binary fluid mixture, where one of the three states 
represents a vacancy instead of a new molecule and the corresponding concentration 
is associated to a variable total volumen. In the last case, the possible phases are 
vapor, miscible liquid mixture and inmiscible liquid mixture. The convenient extensive 
variables are the specific volume v and the relative concentration x = -7- of the two 

c x+y 

fluids and the intensive variables are the pressure and the chemical potential of one of 
the fluids. Moreover, the thermodynamical potential Eq. ([!]) depends on T, v and x 
and is therefore the Helmholtz potential F(T,v,x). This system is perhaps the most 
interesting for applications, given the great amount of experimental data on binary 
fluid mixtures |12|, However, the potential (|IJ) is not the most popular for fitting 
data; a related form which has similar dependence on the relative concentration of the 
two fluids but is of Van der Waals type for the volume is usually considered instead. 
We believe that this form, which is much more difficult to analyze, gives essentially the 
same qualitative behavior. 



4 Applying the CT program 

From Eqs. ([3D and (0), we have a function depending only on two variables, x and y, 
namely, 

H(x, y, a, b, c) = ay(l— x— y)+bx(l— x— y)+cxy+x In x+y\ny+(l — x— y) ln(l— x— y). 

(4) 

Now consider the function H y (x, y, a, b, c) — //2 (where the subindices indicate deriva- 



6 



tives with respect to the variable explicitely written and ix 2 = tx y ) and let the point 
{xqi yo, clq, bo, Co, fx®) be such that H y (P ) — fx® = and H yy (P ) > 0, where P = 
(x ,y ,a ,b Q ,c ). The first condition is the equilibrium condition for y and the sec- 
ond one is required by stability in the y direction. Then, by the Implicit Function 
Theorem, there exists a unique function ip(x, a, b, c, ^2) defined in a neighborhood of 
(xo, ao, bo, Co, fx®) with values in a neighborhood of yo such that H yy {x, y, a, b, c) > 0, in 
these neighborhoods, and H y (x, ip(x, a, b, c, ^2), a, b, c) — /i2 = in the domain of ip and 
4>(xo,a ,bo,co,fxl) = y . 

By solving for y as a function of /i2 and substituting into H—fx 2 y we have performed 
a Legendre transformation, effectively eliminating the variable y. Next we substract 
fx\x (where fx\ = fx x ) to obtain the function 

L(x, a, b, c, jUi, /x 2 ) = # (a;, ^(x, a, 6, c, /i 2 ), a, 6, c) - fx 2 ip(x, a, b, c, /i 2 ) - fx\x, (5) 

representing a one variable non-equilibrium Gibbs potential. In order to have a function 
defined in a neighborhood of = (0, 0, 0, 0, 0, 0) we consider the new function 

L x (x, a, b, c, /ii, fx 2 ) = L(x+x , a+a , b+b , c+c , 1x1+ ix[, fi2+f4)- L ( x o, «o, b , c , fx\, fx®). 

(6) 

Then we have the following equations: 

L 1 (0) = (7) 

and 

L hx (x, a,b, c,fxi,fx 2 ) = 



x + x , ip(x + x , a + a , b + b ,c + c , fx 2 + f4),a + a , b + b ,c + c 



+ H y fa - {fx 2 + H° 2 ) fa - (fix + i4) (8) 

so that 

Li >x {0) = H X {P ) - //?. (9) 

Suppose now that H X (P ) — fx® = (the remaining equilibrium condition). Then 
would be called a critical point of Li(x, 0, 0, 0, 0, 0) (in mathematical terminology). 
This function will be denoted by g\{x), which is the germ to be studied. Of course, 
gi(0) = and g[(0) = from (^) and (|9]). The two first derivatives of g\ are 

g[(x) = H x (x + x , ip(x + x , a , b , c , fx® 2 ), a , b , c ) - //? (10) 

and 

g'((x) = H xx (x + x , fax + x , a , b , c , fx® 2 ), a , b , c ) + H xy (-)fa(-). (11) 
In particular 

g'{(0) = H xx (P ) + H xy (P )iP x (x ,a ,b ,c ,Li 2 v ) (12) 

with 

/ / 7 0\ HxyyPo) /ion 

fa{x , a , b , Co, li 2 ) = - , (13) 
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as deduced from H y — jj 2 = by taking the derivative with respect to x. Suppose that 
the Hessian of H is such that 

H xx (Po) H yy (P ) - H 2 xy (P ) = 0, (14) 

then <7i(0) = and is a degenerate critical point of g±. Finally, we also assume that 
<7i'(0) = ^(O) = (?i(0) = gl l {0) = 0. Then we have imposed five conditions on g\ 
altogether and we should be able to solve for (x , y , a , b , c ). In this case, we say 
that we have reached the highest codimension. We will see in the next section that 
there is indeed such a solution. 



5 Results 



5.1 Highest singularity 

Now we look for a point which fulfills the 5 conditions for the highest singularity 
mentioned above. We succesively have that 



a = c 



2/0 = 1 2 Xq, 

b = Xq 1 , 

1 + 2x 



8x (l-2x ) 



(15) 
(16) 
(17) 



and 



36 4 + 4x0- 1 = 0. 
Thus from the last equation we have that xq = v ^" 1 ~ 0.1201265. Moreover, 

1472 



27 



Xn 



where 



5.C 



256 
162 



' Xr 



157 

To" 



; 2 -17(l-2x ) 2 + 5 3X ° 4X0+1 



x 2 q{1 - 2x )' 



Now we apply results of Singularity Theory || [14 



The 7-jet of g x is j 7 (#i) = ^ gf (0) x 7 with gf(0) ^ 0. 

the essence of g\ respect the identity is 7 and therefore cx(gi) > 7, where o{g\) is 
the determinacy of g\ (see Appendix for the definition of this concept). 

The codimension of j 7 (gi) is cod(j 7 (gi)) = dim vect((x)/(x 6 )) = 5 (see Ap- 
pendix) and o-(j 7 (gi)) < 7. Thus j 7 (gi) is 7-determinate, j 7 (gi) ~ gi and 
gi is 7-determinate. Moreover, cod(#i) = 5 and o"(gi) = 7. A basis of this 
quotient vector space is given by the set {[x],..., [x 5 ]}. Then g±(x, Ai, A5) = 
gi(x) + Xix + X 2 x 2 + ... + A 5 x 5 is a canonical unfolding /c-transversal of the germ 
gi for every k > 0; in particular, for k = 7. Finally g\ ~ x 7 . 
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The L\ function is an unfolding 7-transversal of g\ because one can prove (see 
Appendix) 

(x) = {x 6 )+V Ll + (x) 7+1 , (19) 
where Vl 1 is the real linear space generated by 

{L 1>a {x, 0) - L lia (0, 0), . . . , L 1>fl2 (x, 0) - L^ 2 (0, 0)}, 

where the subindices a, . . . , /x 2 denote derivatives with respect to the correspond- 
ing parameters. Note that Eq. ([H^) has the following expression 

(x) = (x 6 ) + jx (-1 + 2x ) + x 2 -Xq 1 (1 - 3x ) + x 3 (Jxq 1 + ^ 2 ("1 + 3z )) + 

a( 2 8 / x _o-352\ , 1 1 352 , 16 , 

x {- x ° 9 + ^ " X ) x ° 15^ J + ^ (Vfe(3x - 1)^ + jj^o + -x 3 

1 2 ,2 , 4 / 4 _ 2 1 352 _ 2 \ - / 1 352 1 , 

"3 X +X 9^" +X H* +2415-" J +X V24 45- X °~ " X °5!^ 



m o n 2 2 3- 2 -i 4 1 128 _ 2 5 If . 352 
x(l - 2x ) - -x + x — x + x 7^^r x o + x {xq^sx - —x Q 

2 ^ —1 3 —-2 4 —3 5 / / \ 

— X, X -X Q — X -X Q + —X X + —X {—1p5x) 

(20) 

where the Taylor expansions in Vl 1 have been truncated at x 5 because the terms 
of degree sixth and higher are already included in (x 6 ). The resolution of Eq. 
(pOD amounts to prove that a generic nth-degree polynomial with no independent 
term can be generated as a linear combination of the five polynomials between 
the curly brackets. Hence it implies the resolution of a linear system whose 
determinant is x (-58944 xg + 729 ip 5x ) ^ 0. 

Finally, this equality, Eq. (p0|), holds and consequently gi and L\ are isomorphic 
as unfoldings and L\ can be qualitatively studied by the polynomial 

x 7 + Ai x + A2 x 2 + A3 x 3 + A4 x 4 + A5 x 5 , 

which in the terminology of CT corresponds to the Wigwam or Aq catastrophe. 

As a result of this analysis we can affirm that there are three changes of coordinates 
<fi, ip2 and <£>3 and a perturbation e of parameters such that 

Li(x, a, b, c, fix, /i 2 ) = u 7 + Ai« + Xu 2 H h Au 5 + e(a, 6, c, /ii, /i 2 ) (21) 

with 

u = v? 3 pi (fx(x,a,b, 0,111,1^2), 
tp 2 (a, b, c, /ii, /i 2 ) = (Ai, . . . , A 5 ). 
where pi means the first projection, p\ : IR 1+5 — > 1R, that is, Pi(xi, . . . , x 6 ) = x\. 
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Moreover, the bifurcation set of Eq. (|6|) and that corresponding to 

u 7 + X 1 u + Xu 2 + ■■■ + Am 5 (22) 

are diffeomorphic and we rather work with Eq. (p2|) due to its simplicity. The equilib- 
rium manifold in (u, Ax, • • • , A5) is obtained from Eq. ( p2j) by equating the first derivative 
to zero 

7u 6 + Ai + 2A 2 u + 3A 3 u 2 + 4A 4 u 3 + 5A 5 u A = 
and, furthermore, instability occurs if the second derivative also vanishes, 

42u 5 + 2A 2 + 6A 3 u + 12A 4 u 2 + 2OA5 u 3 = 0. 

Now the bifurcation set is obtained by a projection onto the parameter space, that is, 
by eliminating the variable u in this system of two equations. 



5.2 Critical azeotropy as a singularity 

In the process of solving the equations that lead to the highest singularity one goes 
through singularities of lower codimension. They have no particular interest by them- 
selves except in one case, which we proceed to describe. 

Thus we analyze now the five conditions one by one. The equilibrium condi- 
tions allow one to express Xq and yo as functions of the parameters, resulting in 
y = i/)(xq, a , b , c , /i 2 ), already used to define gi(x), and an equation for x derived 
from equating @ to zero. However, we prefer to keep x and y in the equations to 
follow, for it is simpler, understanding that they are to be substituted in the end. 

The first condition on gi(x) is g"(0) = or 

det H i:i = H xx (P ) H yy (P ) - H 2 xy (P ) = 0. (23) 

The Hessian matrix is 

(^)=f" 3 + ai ," 3 V (24) 



with a3 = c — a — b + z 1 , a± = a — c — b + x 1 y ct2 = b — a — c + y 1 . We have 
reinstated z — 1 — x — y for the sake of symmetry and we suppress the subindices zero 



relative to Pq in the next equations. From Eq. ( f24"|) 



det Hij = a 3 «i + 0:3 a 2 + a.\ a 2 = 0, (25) 

which is a quadratic equation on either a, b, c or x, y. The next condition ^"'(O) = is 
equivalent to 

^«2 + ^« 3 --L(a 1 + a 2 ) 3 = (26) 

or 

— a 3 - — {a 3 + ai ) 3 + — al = (27) 
x z y z z A 
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or 

-^(«3 + a 2 ) 3 + ^ + ^ = _ (28) 

Each of these is an equation of third degree in a, b, c. 

Next we consider g™(0) = 0. Again this equation can be separated into three 
symmetric options. Each one imposes new conditions on «j and hence on the elements 
of the Hessian matrix. Furthermore, they can only be fulfilled if two new quantities 
vanish. In other words, the codimension increases by two units with only one condition, 
a non-generic situation. We choose the solution a 3 = a,\ = such that H XX (P ) = 
H xy (P ) = 0. There are two more solutions, obtainable by cyclic permutation of the 
labels of the alphas. As a counterpart of the previous extra increase of the codimension 
by one unit, the next condition is identically fulfilled, gl(0) = 0. Finally, from g]"(0) = 
one obtains the solution quoted in the subsection above. 

Let us compare the foregoing analysis with the one in the previous literature on 
this model |T0| , |2|. This analysis goes as follows. The first condition g"(0) = or, 



equivalently, that the Hessian of H be null has two types of solutions, namely, a simple 
solution, a^ 1 + 1 + ctjT 1 = 0, and a second solution that requires that two of the 
«i vanish simultaneously. The first type is more generic but, unfortunately, does not 
lead to a high codimension singularity, since the equation g l ±(0) = has no solution 
for it. The second solution is the one we have considered in the previous subsection. It 
is called a symmetric solution because the additional condition for the Hessian matrix 
elements implies a symmetry in the phase diagram, namely, a = c. 

Alternatively, we can avoid making any choice on the type of solution until the 
last moment, namely, when we demand g™(0) = 0. Then this condition implies by 
itself (?i(0) = 0. In other words, the singularities given by these two conditions are 
inextricably linked: the first one entails the second one. Moreover, the first condition 
led us to take H XX (P ) = H xy (P ) = in addition to det Hy = 0. We may recall here 
that the vanishing of these two elements of the stability matrix for a binary fluid mixture 
has a thermodynamical interpretation: It occurs when there is critical azeotropy fl2| , 
pages 197-199]. Azeotropy is not a singularity on its own but when it superposes on 
a critical point it enhances its singularity producing a new one. We see what kind of 
singularity it is in our case, namely, the one given by ^[(O) = — a tricritical point. 
At this moment, we are not able to say if this is just a peculiarity of the particular 
three-component model we study or in fact constitutes a general feature. 



6 Discussion 

The solution with the highest codimension, namely, five, which has been found above 
must be isolated; that is, it cannot belong to a continuos family of solutions. Neverthe- 
less, there can be a discrete set of solutions. In fact, we can obtain two other solutions 
from this one by using the symmetry of the original potential G(x,y, z,a,b,c) under 
simultaneous permutations of (x,y,z) and (a, b, c). The way in which they arise is 
clear in the last subsection. It is clear as well that there are no further solutions with 
codimension five. 
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Some readers may be concerned by the fact that our solution is a Landau potential 
which is not bounded below. This does not mean that there is no absolute minimum 
and the Landau potential is meaningless. We must remember that the Landau potential 
is a local object and provides no information on the behavior of the thermodynamical 
potential far from the point that we have called Pq. This point corresponds to a 
degenerate instability rather than to a multicritical point. To be precise, it arises as a 
tricritical point, with Landau potential x 6 , merges with an unstable equilibrium state 
and disappears, for which the appropriate name is tricritical unstable point. Since 
9\ n {fy 7^ 0) there is no tetracritical point in the phase diagram. In this we disagree 
with [jj, where they assert to have three tetracritical points with gl u (0) = 0. We 
attribute this disagreement to a slip on their part, since the conditions they obtain for 
their tetracritical points are precisely the same five conditions we have for our tricritical 
unstable point and the codimension of a generic tetracritical point is six. Thus it seems 
that they are calling tetracritical points to what actually are tricritical unstable points. 

The topology of the overall phase diagram is formed by gluing 3 patches correspond- 
ing each to the phase diagram of the Wigwam catastrophe, as given in the literature 



Tq| . The Taylor series expansions on which the CT is based are supposed to be valid 
in each patch. The coordinates on each patch are not related and therefore we can 
only obtain topological information. Even the topological matching is not a trivial 
matter when we deal with high-dimensional spaces. Since the highest singularities in 
our analysis are essentially the same as those in [[J , so is the phase diagram. The phase 
diagram is of great utility for experimentalists and we should like to obtain some more 
concrete information. We must notice that in a particular model, for example, a bi- 
nary fluid mixture, the energy parameters are fixed and one can only tune the chemical 
potentials — one of which is to be interpreted as pressure in the binary mixture — and 
the temperature. Therefore, one is interested in three-dimensional sections of the over- 
all phase diagram. In these sections it is not generic to have a tricritical point, (which 
has codimension four,) a fact well known to experimentalists. Generic and non generic 
sections of the phase diagram of the Wigwam catastrophe are expected to cover all the 
possibilities. 

Another solution of det Hij = is, of course, that all the matrix elements be null 
so the matrix has rank zero or co-rank 2. This means that one cannot use the Implicit 
Function Theorem to reduce to a function of one variable and one is to proceed with 
the CT program for a function of two variables. The next step is to analyze the 3-jet 
of this function, that is, the form given by the third derivatives. According to its 
signature, given by the sign of its discriminant, there are three possible cases: If it is 
negative, the canonical form is x 3 — 3xy 2 , called the elliptic umbilic catastrophe; if it 
is positive the canonical form is x 3 + 3 x y 2 , called the hyperbolic umbilic catastrophe; 
if it is null, the 3-form is degenerate, and one is to anlyze the 4-jet to determine the 
type of singularity, which may be the parabolic umbilic catastrophe or a type even more 
complex. The calculations driving at establishing the highest singularity in our five 
dimensional parameter space are complicated and the results for co-rank two shall be 
reported in the future |L7|. 
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A Summary of Catastrophe Theory 



In this Appendix we are going to present the main mathematical concepts widely 



introduced in Refs. |JL4| , |15| and necessary to follow the main steps developed in Section 
5. 

Let us consider real functions of class oo and defined in a neighbourhood of G IR n . 
We establish that two functions are equivalent if they coincide in a neighbourhood 
of 0. The classes we obtain are called germs of functions and the set of germs is 
denoted by E(n). The operations / + g and / • g give to E(n) the estructure of a 
ring and M(n) — {./ G E(n)/f(0) = 0} is a maximal ideal of this ring. Moreover, 
the operations f + g and A ■ / with A G IR give to E{n) the structure of a real 
vector space of dimension oo. The ideal M(n) is generated by x±,...,x n , that is, 

M ( n ) = {fx%i H h fnX n I fx, ■ ■ ■ , fn e E(n)}. In general, if ft, ...,/„ G E(n), we 

designate by (fi, . . . , / n ) to the ideal generated by fi, . . . , /„, that is, 

(fx,---,fn) = {fl9l + ■■■ + fn9n / 9x, ■ ■ ■ , 9n £ 

In particular, M(n) = (x±, . . . , x n ). 

It is possible to define powers of M(n) as M(n) k . It can be proven that M(n) k 
is equal to the ideal of E{n) generated by the monomials in x\, . . . ,x n of degree k. 
In particular, for example, (x,y) 2 = (x 2 ,xy,y 2 ). We have also that M(n) k+1 = {/ G 
E(n)/D 1 f(0) =0,i < k} where with D l we mean the derivative of degree i. 

In a similar way, we can define the germs of dipheomorphisms of class oo from ]R n 
in R" which transform in 0. This set is denoted by G(n). We say that two germs 
are equivalent if there exists a change of coordinates (f G G(n) such that / = g if and 
it is denoted by / ~ g. When the fc-jets (k G IN )of two functions are equal we say 
that these functions are /^-equivalent (/ ~& g). A germ f is fc-determinate if for every 
germ g such that both fc-jets are equal we have that / ~ g. The determinacy of a germ 
f is the smallest number k G IN such that / is /c-determinate and it is called cr(/). 
Therefore, we notice that : 

If f is fc-determinate then / ~ j k f , 



If f is fc-determinate and / ~^ g then g is /c-determinate , 
If f is fc-determinate and / ~ g then g is fc-determinate. 

The ideal of Jacobi of a germ f is defined by 

A(f) = (D x J,...,D x J) 

where D X1 , . . . , D Xn are the partial derivatives with respect to the xi, . . . , x n variables. 
If/~<7then A(f) = A(g). 

Now let us assume that / G M(n) 2 . Then A(/) C M(n) and we can speak of 
the quotient vector space M(n)/A(/). The dimension of this vector space is called 
codimension of /, cod(/). It can be proven that if / G M{n) then cod(/) is finite 
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if and only if a(f) is finite and in this case a(f) — 2 < cod(/). Moreover, if / ~ g 
then cod(/) = cod(g). Thus the codimension generally coincides with the number of 
conditions necessary to specify a function germ up to diffeomorphisms, which is the 
usual geometrical concept. 

If g G M(n + r) and / G M(n) 2 , we say that g is a r-unfolding of / if g(x, 0) = f(x) 
(with r parameters which we define as yi, . . . , y r ). Now let us assume that g G M(n+r) 
is an unfolding of / G M(n) 2 and k G INT. We say that g is A;-transversal if 

M(n) = A(f) + M(n) k+1 + V g 

where V g is the real vector subset of M{n) generated by the vectors D yi g(x,0) — 
A/ifi^O, 0), . . . , D Vr g(x, 0) — Dy r g(0, 0). 

Finally, the theorem of A;-transversality for unfoldings can be stated as follows: 
Let us consider / G M(n) 2 /c-determinate and g and h two unfoldings of / with r 
parameters which are fc-tranversal. Then g and h are isomorphic. 
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